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Using the OPAL detector at LEP, the running of the effective QED coupling a(t) is measured 
for space-like momentum transfer, 2 < — t < 6 GeV 2 , from the angular distribution of small- 
angle Bhabha scattering. This is currently the most significant direct observation of the 
running of the QED coupling in a single experiment and the first clear evidence of the hadronic 
. contribution to the running in the space-like region. Our result is in good agreement with 

f*^ ' standard evaluations of a(t), based on data in the time-like region. 

in 
o 
in 
o 



> 



- 1—1 

X 



Introduction 



The effective QED coupling a(t) is an essential ingredient for many precision physics predictions. 
It contributes one of the dominant uncertainties in the electroweak fits constraining the Higgs 



mass. The effective QED coupling is generally expressed as: 



Q < 4 > = T^tf) « 



mm® 



where Qo = a(t = 0) — 1/137 is the fine structure constant, t is the momentum transfer 
squared of the exchanged photon and Aa is the vacuum polarization contribution. Whereas the 
leptonic contributions to Aa are calculable to very high accuracy, the hadronic ones have to 
be evaluated by using a dispersion integ ral o ver the measured cross section of e + e~ — > hadrons 
at low energies, plus perturbative QCD^^. There are also many evaluations which are more 
theory-driven, extending the application ofperturbative QCD down to ~ 2 GeV (see for example 
the reference'"'). An alternative approach^ uses perturbative QCD in the negative t (space-like) 
region. 

There have been only a few direct observations of the running of the QED coupling'-' 
Here we present a new result from the OPAL collaboration. A full description can be found 
in the OPAL paper ®. The running of a is measured in the space-like region, by studying 
the angular dependence of small-angle Bhabha scattering, e + e~ — * e + e~, at LEP. Small-angle 
Bhabha scattering appears to be an ideal process for a direct measurement of the running of 
a(t) in a single experiment, as it is an almost pure QED process, strongly dominated by t- 
channel photon exchange. Moreover the data sample has large statistics and excellent purity. 
The Bhabha differential cross section can be written in the following form for small scattering 
angle: 

da d<r(°) fa(t)\ 2 M . M . , . 



where da^ /dt = Afta^/i 2 is the Born term for the t-channel diagram, e represents the radiative 
corrections to the Born cross section, while <5 7 and 5z are the interference contributions with s- 
channel photon and Z exchange respectively. <5 7 and 5% are much smaller than e and the vacuum 
polarization. Therefore, with a precise knowledge of the radiative corrections (e term) one can 
determine the effective coupling a(t) by measuring the differential cross section. This method 
has also been advocated in a recent paper^l 



2 Detector and event selection 



We use OPAL data collected in 1993-95 at energies cl ose to the Z resonance peak. In particular 
this analysis is based on the OPAL SiW luminomet er 1*U The SiW consisted of two cylindrical 
calorimeters encircling the beam pipe at a distance z ~ ±2.5 m from the interaction point. Each 
calorimeter was a stack of 19 silicon layers interleaved with 18 tungsten plates, with a sensitive 
depth of 14 cm, representing 22 radiation lengths (Xo). The sensitive area fully covered radii 
between 6.2 and 14.2 cm from the beam axis, corresponding to scattering angles between 25 and 
58 mrad. Each detector layer was segmented with R-<f> geometry in a 32 x 32 pad array. The 
pad size was 2.5 mm radially and 11.25 degrees in azimuth. In total the whole luminometer 
had 38,912 readout channels corresponding to the individual silicon pads. Particles coming from 
the interaction point had to traverse the material constituting the beam pipe and its support 
structures as well as detector cables before reaching the face of the SiW calorimeters. This 
preshowering material was minimum near the inner angular limit, about 0.25 Xo, while in the 
middle of the acceptance it increased to about 2X(j. When LEP2 data-taking started in 1996 
the detector configuration changed, with the installation of tungsten shields designed to protect 
the inner tracking detectors from synchrotron radiation. This reduced the useful acceptance of 
the detector at the lower angular limit. Therefore we limited this analysis to the LEP 1 data. 

The event selection is similar to the one used for luminosity measurements^^. The selected 
sample is strongly dominated by two-cluster configurations, with almost full energy back-to-back 
e + and e _ incident on the two calorimeters. At leading order the momentum transfer squared t 
is simply related to the scattering angle 6, which is measured from the radial position R of the 
scattered e + and e~ at reference planes located within the SiW luminometers: 

1 - cos 6 s6 2 . , . 

t = — s — — ; tan = R/z . (3) 

At the center-of-mass energy y/s ~ 91 GeV our angular acceptance corresponds to 2 < — t < 
6 GeV 2 . 

The radial distributions are shown in Fig.^for the complete data statistics, compared to the 
Monte Carlo distributions normalized to the same number of events. Due to the back-to-back 
nature of Bhabha events, the two sides do not contribute independent statistical information. 
After the studies mentioned in section 0J we decided to use the Right side distribution for the 
final fits, to keep at minimum possible unassessed systematic errors. Consistent results are 
obtained with the use of the Left side distribution. 



3 Fit method 

The counting rate of Bhabha events in the SiW is used to determine the integrated luminosity, 
so that we cannot make an absolute measurement of a(t) without an independent determination 
of the luminosity. 

We compare the radial distribution of the data (and hence the t-spectrum) with the predic- 
tions of the BHLUMI Monte Carlo ±z . This is a multiphoton exponentiated generator accurate 
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Figure 1: Radial distributions for the complete data statistics. The points show the data and the histogram the 
Monte Carlo prediction, assuming the expected running of a, normalized to the same number of events. The 

lower plots show the ratio between data and Monte Carlo. 



up to the leading logarithmic 0(a 2 L 2 ) terms . Higher order photonic contributions are partially 
included by virtue of the exponentiation. It has been used to determine the luminosity at LEP 
and has been widely cross-checked with many alternative calculations. If the Monte Carlo is 
modified by setting the coupling to the constant value a(t) = «o, the ratio / of the number of 
data to Monte Carlo events in a given radial bin is: 

/(*) = ArO 7* K 1 ^77T • ( 4 ) 



JV& (t) U-Aa(t) 

The dominant dependence of Aa(t) expected from theory is logarithmic. We therefore fitted 
the ratio f(t) as: 

f(t) = a + bln(f) (5) 



where to = —3.3 GeV 2 is the mean value of t in the data sample. The parameter a, about 
unity, is not relevant since the Monte Carlo is normalized to the data. The slope b represents 
the full observable effect of the running of a(t), both the leptonic and hadronic components. It 
is related to the variation of the coupling by: 

Aa(t 2 ) - Aa(ti) ~ - In (-) (6) 



2 \*i 

where t\ = —1.81 GeV 2 and t% = —6.07 GeV 2 correspond to the acceptance limits. 
4 Main systematic effects 

It is important to realize which systematic effects could mimic the expected running or disturb 
the measurement. The most potentially harmful effects are biases in the reconstructed radial 

a L = ln(|t|/m.e) — 1 is the large logarithm. 



Table 1: Fit result for each dataset and average. For each value of b the first error is statistical and the second 

the full experimental systematic. 
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coordinate. Most simply one could think of dividing the detector acceptance into two and 
determining the slope using only two bins. In such a model the running is equivalent to a bias in 
the central division of 70 /im. Biases on the inner or outer radial cut have a little less importance 
and could mimic the full running for 90 or 210 fim systematic offsets respectively. Concerning 
radial metrology, a uniform bias of 0.5 mm on all radii would give the same observable slope 
as the expected running. Knowledge of the beam parameters, particularly the transverse offset 
and the beam divergence, is also quite important. Thus, limitation of systematic error in the 
reconstructed radial coordinate is key to the current measurement. 

Details of how the coordinates are formed from the recorded pad information are found in 
The fine radial and longitudinal granularity of the detector are exploited to produce precise 
radial coordinates. The reconstruction determines the radial coordinate of the highest energy 
cluster, in each of the Right and Left calorimeters. Each coordinate uses a large number of 
pads throughout the detector, from many silicon layers, and is projected onto a reference layer, 
close to the average longitudinal shower maximum. The residual bias, or anchor, of this radial 
coordinate is then estimated at each pad boundary in a given layer of the detector. Here we 
rely on the fact that, on average, the pad with the maximum signal in any particular layer will 
contain the shower axis. Then from the anchors we obtain bin-by-bin acceptance corrections 
which are applied to the radial distribution. This procedure, named anchoring, is the most 
delicate part of the analysis, and was carefully studied®. The challenging aspect is controlling 
the residual bias on the radial coordinate to a level below ~ 10 /im uniformly throughout the 
acceptance. 



5 Results 

The ratio of data to Monte Carlo is fitted to Eq. and the results are given in Table The 
nine datasets give consistent results, with x 2 /d.o.f = 6.9/8 for the average b considering only 
statistical errors. The most important systematic errors come from the anchoring procedure 
and the preshowering material, both affecting the radial coordinate. The fit results are then 
combined, by considering the full error correlation matrix, obtaining: 



b = (726 ± 96 ± 70 ± 50) x 10~ 5 
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Figure 2: |t| spectrum normalized to the BHLUMI theoretical prediction for a fixed coupling (Aa = 0). The 
points show the combined OPAL data with statistical error bars. The solid line is our fit. The horizontal line 
(Ratio=l) is the prediction if a were fixed. The dot-dashed curve is the prediction of running a determined by 
vacuum polarization with only virtual e + e~ pairs, the dashed curve includes all charged lepton pairs and the 
dotted curve the full Standard Model prediction, with both lepton and quark pairs. 

where here, and also in the results quoted below, the first error is statistical, the second is the 
experimental systematic and the third is the theoretical uncertainty. The total significance of 
the measurement is 5.6 a. 

The theoretical uncertainty is dominated by the photonic corrections to the leading i-channel 
diagram, in particular by missing 0(a 2 L) terms, and the technical precision of the calculation. 
We estimated these uncertainties by comparing BHLUMI with alternative Monte Carlo calcu- 
lations. Other uncertainties, from Z interference and the contribution of light e + e~ pairs, were 
also estimated and added in quadrature. 

The result for the combined data sample is illustrated in Fig. El The logarithmic fit to 
Eq. [5] describes the data very well, x 2 /d.o.f = 1.9/3, although a simple linear fit would also be 
adequate, giving x 2 /d.o.f = 2.7/3. The data are clearly incompatible with the hypothesis of 
a fixed coupling. The fitted logarithmic dependence agrees well with the full Standard Model 
prediction including both leptonic and hadronic contributions, with the hadronic part obtained 

lol 

by the Burkhardt-Pietrzyk parameterization 1 ™ 

The effective slope gives a measurement of the variation of the coupling a(t) from Eq. |BJ 

Aa(-6.07 GeV 2 )- Aa(- 1.81 GeV 2 ) = (440 ± 58 ± 43 ± 30) x 10~ 5 . 

This is in good agreement with the Standard Model prediction, which gives 5 (Aa) = (460 ±16) x 
10 -5 for the same t interval, where the error originates from the uncertainty of the hadronic 




component. The evaluation^ of Acth nri h as a relative precision ranging from 2.5% at t = 
-1.81 GeV 2 to 2.7% at t = -6.07 GeV 2 Ql. 

The absolute value of Aa in our range of t is expected to be dominated by e + e~ pairs, with 
the relevant fermion species contributing in the approximate proportions: e : fx : hadron ~ 4 : 
1 : 2. Our measurement is sensitive, however, not to the absolute value of Aa, but only to its 
slope within our t range. Contributions to the slope b in this range are predicted to be in the 
proportion: e : fi : hadron ~ 1 : 1 : 2.5. Fig. [21 shows these expectations graphically. We can 
discard the hypothesis of running due only to virtual e + e _ pairs with a significance of 4.4 a. 

The data are also incompatible with the hypothesis of running due only to leptons. If we 
subtract the precisely calculable theoretical prediction for all leptonic contributions, 5(Aa\ ep ) = 
202 x 10 -5 , from the measured result, we can determine the hadronic contribution as: 

Aa had (-6.07GeV 2 ) - Aa had (-1.81 GeV 2 ) = (237 ± 58 ± 43 ± 30) x 10~ 5 . 

This has a significance of 3.0 a, considering all the errors. 

Our result can be easily compared to the previous one by L3^. If the latter is expressed as 
a slope according to Eq. El it becomes: b^ L3 ^ = (1044 ± 348) x 10 -5 . The two measurements are 
shown in Fig. |31 The L3 result has a larger error dominated by experimental systematics but is 
consistent with ours. The average gives: iA ave ) = (759 db 113 db 50) x 10~ 5 , where the first error 
is obtained from the experimental errors and the second is the theoretical uncertainty that we 
estimated for our measurement, which will likely be common. The average is in good agreement 
with the prediction using the Burkhardt-Pietrzyk parameterization. 
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Figure 3: Effective slope 6 = 2 < dAa/dlnt >. The OPAL and L3 measurements are shown together with their 
average. The solid line is the SM prediction, with the band showing its uncertainty. The dashed line at b = 

represents the case of no running. 



6 Conclusions 

We have measured the scale dependence of the effective QED coupling from the angular dis- 
tribution of small-angle Bhabha scattering at LEP, using the precise OPAL Silicon- Tungsten 
luminometer. We obtain the strongest direct evidence for the running of the QED coupling 
ever achieved in a single experiment, with a significance above 5 a. Moreover we report the 
first clear experimental evidence for the hadronic contribution to the running in the space-like 
region, with a significance of 3 a. This measurement is one of only a very few experimental tests 
of the running of a(t) in the space-like region, where Aa has a smooth behaviour. Our result is 
in good agreement with standard evaluations of a(t), based on data in the time-like region. 
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